Possible experimental manifestations of the many-body localization 
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Recently, it was predicted that if all one-electron states in a non-interacting disordered system are 
localized, the interaction between electrons in the absence of coupling to phonons leads to a finite- 
temperature metal-insulator transition. Here, we show that even in the presence of a weak coupling 
to phonons the transition manifests itself (i) in the nonlinear conduction, leading to a bistable I-V 
curve, (ii) by a dramatic enhancement of the nonequilibrium current noise near the transition. 



Introduction. — Low-temperature charge transport in 
disordered conductors is governed by the interplay be- 
tween elastic scattering of electrons off static disorder 
(impurities) and inelastic scattering (electron-electron, 
electron-phonon, etc.). At low dimensions an arbitrarily 
weak disorder localizes Q all single-electron states 0, HI , 
and there would be no transport without inelastic pro- 
cesses. For the electron-phonon scattering, the dc con- 
ductivity o~(T) at low temperatures T is known since long 
ago in d dimensions 



lncr(T) cx -1/T 7 , 7 = l/(d + l). 
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What happens if the only possible inelastic process is 
electron-electron scattering? The answer to this ques- 
tion was found only recently 0: a(T) = identically for 
T < T c , the temperature of a metal-insulator transition. 
Here we discuss experimental manifestations of this tran- 
sition in real systems, where both electron-electron and 
electron-phonon interactions are present. We show that 

(i) the I-V characteristic exhibits a bistable region, and 

(ii) non-equilibrium current noise is enhanced near T c . 
The notion of localization was originally introduced 

for a single quantum particle in a random potential 
Subsequently, the concept of Anderson localization was 
shown to manifest itself in a broad variety of phenomena 
in quantum physics. This concept can also be extended to 
many-particle systems. Statistical physics of many-body 
systems is based on the microcanonical distribution, i. e., 
all states with a given energy are assumed to be realized 
with equal probabilities. This assumption means delocal- 
ization in the space of possible states of the system. It 
does not hold for non-interacting particles; however, it is 
commonly believed that an arbitrarily weak interaction 
between the particles eventually equilibrates the system 
and establishes the microcanonical distribution. 

Many-body dynamics of interacting systems and its re- 
lation to Anderson localization has been discussed in the 
context of nuclear Q and molecular [3] physics. For in- 
teracting electrons in a chaotic quantum dot this issue 
was raised in Ref. [U, where it was shown that electron- 
electron interaction may not be able to equilibrate the 
system. This corresponds to Anderson localization in the 
many-body space. Recently it was demonstrated that 
in an infinite low-dimensional system of (weakly) inter- 



FIG. 1: Schematic temperature dependence of the dc con- 
ductivity cr(T) for electrons subject to a disorder potential 
localizing all single particle eigenstates, in the presence of 
weak short-range electron-electron \S^, A < 1, established in 
Ref. pj. Below the point of the many-body metal-insulator 
transition, T < T c ~ (5^/ [ A In A| , no inelastic relaxation occurs 
and cr(T) =0. At T ^> <5<;/A 2 the high-temperature metallic 
perturbation theory^] is valid, and corrections to the Drude 
conductivity <7oo are small. In the interval 8^/\ <C T <C S^/\ 2 
electron-electron interaction leads to electron transitions be- 
tween localized states, and the conductivity depends on tem- 
perature as a power-law. 



acting electrons, subject to a static disorder, Anderson 
transition in the many-body space manifests itself as a 
finite-temperature metal-insulator transition [f|. 

Let single-particle eigenstates be localized on a spatial 
scale C/oc (localization length). The characteristic energy 
scale of the problem is the level spacing within the lo- 
calization volume: 8^ = l/(^Czo C )' v being the density of 
states per unit volume. Below we neglect energy depen- 
dences of v, C/oc, and 5q. According to Ref. [5(, as long 
as electrons are not coupled to any external bath (such 
as phonons), a weak short-range electron-electron inter- 
action (of typical magnitude X5^, with the dimensionless 
coupling constant A < 1) does not cause inelastic relax- 
ation unless the temperature T exceeds a critical value: 
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The small denominator | A In A| represents the character- 
istic matrix element of the creation of an electron-hole 
pair. The ratio T c /8q is the number of states available 
for such a pair (in other words, the phase volume) at 
T = T c . Only provided that this number is large enough 
to compensate the smallness of the matrix element, the 
interaction delocalizes the many-body states and thus 
leads to an irreversible dynamics. As a result, the finite- 
temperature dc conductivity <j(T) vanishes identically if 
T < T c , while a(T > T c ) is finite, i. e. at T = T c a metal- 
insulator transition occurs. The overall dependence cr(T) 
is summarized in Fig. [1] 

In any real system the electron-phonon interaction is 
always finite. This makes u(T) finite even at T < T c : 
<t(T) is either exponentially small [Eq. Jl}] at T <C <5<j, 
or follows a power- law [lTj at 5 ( < T < T c . At the 
transition point the phonon-induced conductivity is not 
exponentially small, i. e. phonons smear the transition 
into a crossover. Are there any experimental signatures 
of the many-body localization? In what follows we show 
that if the electron-phonon coupling is weak enough, a 
qualitative signature of the metal-insulator transition can 
be identified in the nonlinear conduction. Namely, in a 
certain interval of applied electric fields £ and phonon 
temperatures T p h both metallic and insulating states of 
the system turn out to be stable! As a result, the I- 
V curve exhibits an S-shaped bistable region (Fig. [3]). 
Moreover, we show that the many-body character of 
the electron conduction dramatically modifies the non- 
equilibrium noise near the transition [Eq. ([12|>] . 

Bistable TV curve. — Our arguments are based on two 
observations. First, in the absence of phonons a weak but 
finite electric field cannot destroy the insulating state - it 
rather shifts the transition temperature. Let us neglect 
the effect of the field on the single-particle wave functions, 
representing it as a tilt of the local chemical potential of 
electrons. Then at T = the role of the field in the insu- 
lating regime is increase the energy of the electron-hole 
(e-h) excitation of a size L by e£L. This provides in addi- 
tional phase volume of the order oieSL/Sc,. However, for 
L > C/ oc the matrix element for creation of such an ex- 
citation quickly vanishes. In the diagrammatic language 
for the effective model of Ref. [H this means that each 
electron-electron interaction vertex must be accompanied 
by tunneling vertices which describe coupling between lo- 
calization volumes and whose number is (i) at least one 
in order to gain phase volume (in contrast to the finite-T 
case when tunneling had to be included only to overcome 
the finiteness of the phase space in a single grain [f| ) , and 
(ii) not much greater than one, otherwise the diagram is 
exponentially small. As a result, at T = the insulator 
state is stable provided that £ < £ c ~ T c /(eQ oc ). 

In the same way one can analyze the finite-temperature 
correction to the critical field, and the finite-field correc- 
tion to the critical temperature can be found by taking 
into account the extra phase volume in the calculation of 



i (T d - T ph )/eS, 




Tph T c (£) T e i 

FIG. 2: Sketch of the dependences (T e! -T ph )/(e£) (dashed 
line) and L p h(T e i) (solid lines, electron-phonon coupling 
strength being weaker for higher curves) . Actual value of T B i 
is determined by their crossing. 

Ref. @. One obtains a(T) = for T <T C (£), where 

T c (£) = T C - cie£Q oc . (3) 

with a model-dependent factor c\ ~ 1, weakly dependent 
on £ [here and below T c without the argument £ is the 
zero-field value given by Eq. (J2J], As a consequence, at 
T < T c {£) the nonlinear transport, as well as the linear 
one, has to be phonon-assisted. 

The second observation is that when both a and £ 
are finite, there is Joule heating. The thermal balance is 
qualitatively different in the insulating and the metallic 
phases. Deep in the insulating phase (T <C T c ) each 
electron transition is accompanied by a phonon emis- 
sion/absorption, i. e. electrons are always in equilibrium 
with phonons whose temperature T p h <C T c we assume 
to be fixed. On the contrary, in the metallic phase elec- 
trons gain energy when drifting in the electric field, i. e. 
they are heated. Due to this Joule heating the effective 
electron temperature T e i deviates from the bath temper- 
ature. The role of phonons is then to stabilize T e i- For 
weak electron-phonon coupling T e i and T p h can differ sig- 
nificantly. A self-consistent estimate for T e i follows from 

T e i - T ph ~ e£L ph (T e i) , (4) 
L ph (T e i ) = \J D(T e i) T p h(T e i) . (5) 

Here T p h(T e i) is the time it takes an electron to emit 
or absorb a phonon, L p h(T e i) is the typical electron dis- 
placement during this time, and D(T e i) = a(T e i)/(e 2 i/) 
is the electron diffusion coefficient. 

We sketch in Fig. (T e i - T ph )/(e£) and L ph {T e i) 
for different electron-phonon coupling strengths as func- 
tions of T e i. It is taken into account that (i) L p h coin- 
cides with variable range hopping length at T e i -C §(;, 
(ii) L ph ~ Cioe at 6 C « Tei « T c (£), (iii) D(T el ) 
quickly rises to its large metallic value Dao ~ ^cC/ 2 oc near 
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FIG. 3: (a) sketch of the bistable I-V curve for a fixed value 
of T p h\ (b) (£, T'ph)-plane with the bistable region schemat- 
ically shown by shading. The dashed line represents the 
crossover between the metallic state at high electric field £ 
or high phonon temperature T p h, and the insulating state at 
low £ and low T p h- 



T c (£), (iv) T p h(T e [) decreases as a power law with in- 
creasing T e i. The peak of the curve rises with decreas- 
ing electron-phonon coupling strength, and eventually 
the curve crosses the straight line. After that, in ad- 
dition to T e \ = Tph, Eq. ([4| acquires two more solutions, 
both with T e i > T c (£), of which only the rightmost so- 
lution is stable. The maximum of L p h can be estimated 
as L* h ~ ^/D 00 t p ^(T c ), so three solutions appear when 
T c (£ ) — Tph <C e£L* h . At the same time, as seen from 
Eq. (|3|), electric field is unable to break down the insu- 
lator as long as e£Q oc <C T c — T p h- Thus, the interval 
of electric fields where both regimes are stable, is deter- 
mined by 



T c (£) - Tph c T ? ~ T P h 
eL ph eQoc 

The two conditions are compatible provided that 

Lph ^ Cl0C 7 



(0) 



(7) 



which is realistic when electron-phonon coupling is weak. 

In the bistable region (J6j) , for a given value of £ one 
finds two stable solutions for T e i, giving two possible 
values of the conductivity and the current, which corre- 
sponds to an S'-shape current- voltage characteristic ll| , 



the third (unstable) solution corresponding to the neg- 
ative differential conductivity branch. The macroscopic 
consequences of such behavior depend on the dimension- 
ality. In a 2d sample the two phases of different electronic 
temperature and current density can coexist, separated 
by a boundary of the width ~ L* h , parallel to the di- 
rection of the electric field. The particular state of the 
system determined by the boundary conditions (proper- 
ties of the external circuit), as well as by the history. 

Noise enhancement. — In the vicinity of the criti- 
cal point conduction is dominated by correlated many- 
electron transitions (electronic cascades). Each cascade 
is triggered by a single phonon. As T e i — ► T c {£), the typ- 
ical value n of the number n of electrons in the cascade 
diverges together with the time duration of a cascade. 
The results of Ref. Q , adapted for a finite electric field, 
give the following probability for an n-electron transition 
to go with the rate T: 



4tt r 3 / 2 



r n oc 



T P , 



cie£(. 



loc 
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which gives 



1/n ~ \n[T c /(T e i + cie£(i oc )]. 



(8) 



(9) 



The divergence in n is cut off when electron-phonon cou- 
pling is finite. The largest n is such that the phonon 
broadening of the single-electron levels, l/T p h(T c ), is 
comparable to fi-particle level spacing (in other words, 
time duration of a cascade cannot exceed r p h ) ■ 
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ph) 



max 4 ln(T c /(J c ) ' 
(10) 

with logarithmic precision; n a represents the di verg ent 
spatial extent of the cascade (correlation length) [12| . 

Each many-electron transition can be characterized, 
besides its rate T, by the total dipole moment d it 
produces. The corresponding backward transition pro- 
duces the dipole moment — d and goes with the rate 
re -(£-d)/T eI ^ The average current (j( t )} is d eter . 

mined by the difference between forward and backward 
rates; obviously, it vanishes for £ — 0. At the same time, 
the noise power (second cumulant) S2 = J[(I(t)I(t')) — 
(I(t))(I(t'))] dt' is determined by the sum of the forward 
and backward rates; at £ = it is given by the equilib- 
rium Nyquist-Johnson expression. 

Equilibrium noise carries no information about the na- 
ture of conduction. To see a signature of many-electron 
transitions it would be natural to analyze the shot noise, 
whose power is proportional to the charge transferred in 
a single event. Many-electron cascades would then cor- 
respond to "bunching" of electrons, thus increasing the 
shot noise. However, shot noise is observed in the limit 
when transitions transferring charge only in one direc- 
tion (namely, d ■ £ > 0) are allowed, i. e. T e i <C e£Q oc , 
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which is impossible to satisfy in the insulating state, 
as T e i ~ max{ T p h, e£Q oc } [14J. Thus, S2 inevitably 
has both equilibrium and non-equilibrium contributions, 
which are difficult to separate. 

To see the "bunching" effect unmasked by a large ther- 
mal noise at low fields one should study the third Fano 



factor S3 of the current fluctuations [15j. Indeed, being 
proportional to an odd power of the current, it vanishes 
in equilibrium, so it is not subject to the problems de- 
scribed above for 52. In a wire of length L the ratio 
S3 /(I) is given by 



5s_ 

(I) 



■ «r*» 



l- 1 «rd» 



(ii) 



The double angular brackets on the right-hand side mean 
the sum over all allowed transitions. For nearest-neighbor 
single-electron transitions with d = ±e(i oc Eq. (TlTj) gives 
53/(1) ~ e 2 ((i oc /L) 2 , which is analogous to the Schot- 
tky expression reduced by the effective number of tunnel 
junctions in series, L/Q oc , for 52 

Since d 3 diverges stronger than d as n — + 00, we ex- 
pect a divergence in Eq. (fTTj) . The critical index of d 
depends on the order of limits: d ~ y/neQ oc if the lin- 
ear response limit £ — > is taken prior to n — > 00, while 
d ~ ne(i oc (e£Qoc/T c ) for a small but finite £. As a result, 



S3 
(I) 



max < n , 



ne£Q, 



(12) 

where n is given by Eq. @, and the saturation of the 
divergence is determined by Eq. (jlOp . Upon further in- 
crease of the temperature, the system crosses over to the 
metallic state, and n starts to decrease. This decrease is 
governed by the same Eq. (fl~0|) with the phonon inelas- 
tic rate substituted by the typical value of the electron- 
electron inelastic rate, which grows with temperature. 
As the critical behavior of T on the metallic side of the 
transition is unknown, we cannot give any quantitative 
estimate of 53 above T c . 

Conclusions. — In conclusion, we have shown that the 
finite-temperature metal-insulator transition, predicted 
theoretically in Ref . [B} , can manifest itself on the macro- 
scopic level as an S-shape current-voltage characteristic 
with a bistable region. In fact, the hysteretic behaviour 
of the current in Y^Sii-a; 17J is a possible candidate for 



the effect discussed in the present paper. 

Besides, we have shown that the many-body nature of 
the conduction near the transition manifests itself in the 
dramatic increase of the non-equilibrium current noise: 
the noise depends on the total charge transferred in each 
random event, while the number of electrons, involved in 
such an event, increases as one approaches the transition. 
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